Event-by-event shape and flow fluctuations of relativistic heavy-ion collision fireballs 
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Heavy-ion collisions create deformed quark-gluon plasma (QGP) fireballs which explode anisotrop- 
ically. The viscosity of the fireball matter determines its ability to convert the initial spatial defor- 
mation into momentum anisotropies that can be measured in the final hadron spectra. A quantita- 
tively precise empirical extraction of the QGP viscosity thus requires a good understanding of the 
initial fireball deformation. This deformation fluctuates from event to event, and so does the finally 
observed momentum anisotropy. We present a harmonic decomposition of the initial fluctuations 
in shape and orientation of the fireball and perform event-by-event (2+l)-dimensional ideal fluid 
dynamical simulations to extract the resulting fluctuations in the magnitude and direction of the 
corresponding harmonic components of the final anisotropic flow at midrapidity. The final harmonic 
flow coefficients are found to depend non-linear ly on the initial harmonic eccentricity coefficients. 
We show that, on average, initial density fluctuations suppress the buildup of elliptic flow relative to 
what one obtains from a smooth initial profile of the same eccentricity, and discuss implications for 
the phenomenological extraction of the QGP shear viscosity from experimental elliptic flow data. 

PACS numbers: 25.75.-q, 12.38. Mh, 25.75.Ld, 24.10.Nz 
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I. INTRODUCTION 



In ultrarelativistic heavy-ion collision experiments, a 
fraction of the incoming kinetic energy is converted into 
new matter deposited in the collision zone. The distribu- 
tion of this matter in the plane transverse to the colliding 
beams is inhomogeneous and fluctuates from collision to 
collision. At the collision energies available at the Rela- 
tivistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC), the produced matter is sufficiently dense 
and strongly interacting that it quickly reaches a state of 
approximate local thermal equilibrium. Its subsequent 
evolution can thus be described by fluid dynamics until it 
eventually becomes too dilute and breaks apart. Hydro- 
dynamic forces (i.e. pressure gradients) convert the inho- 
mogencities and deformations of its initial spatial density 
distribution into anisotropies of the final hydrodynamic 
flow. The latter can be extracted from the momentum 
distributions of the finally emitted particles. The effi- 
ciency with which the geometric deformation and fluctu- 
ating inhomogencitics in the initial density distribution 
are converted into final flow anisotropies is controlled by 
the viscosity of the expanding fluid. For a given source 
deformation, ideal fluid dynamics generates the largest 
flow anisotropy; it corresponds to the limit of zero mean 
free path and instantaneous thermalization, which allows 
for the largest possible collective response, via final state 
interactions, to irregularities in the geometric structure 
of the fireball [1]. Viscosity accounts for finite interac- 
tion cross sections and non-zero mean free paths which 
reduce the amount of flow anisotropy that can be gen- 
erated from a given geometric deformation. In hydrody- 
namic language, viscous pressure components inhibit the 
development of flow anisotropies and tend to smoothen 
irregularities in the flow distribution. By measuring the 
flow anisotropies and relating them to the initial geomet- 
ric deformations (as calculated from theoretical models 



for the collision geometry) one can, in principle, deter- 
mine the fluid's viscosity experimentally [2-6]. 

Until recently, most of the attention has been focussed 
on elliptic flow i> 2 = (cos(20 p )) and its relation to the spa- 
tial eccentricity e x = j^r^sj = (r 2 cos(2</> 3 ))/ (r 2 ) (more 

precise definitions will be given in Sec. II) [2-12]. Event- 
by-event fluctuations in the initial state were only treated 
on average, by taking into account their effects on the av- 
erage eccentricity of an ensemble of events [5, 6, 11-13] 
and then propagating a smooth initial density profile cor- 
responding to that ensemble average hydrodynamically. 
In this way one can only compute the average elliptic 
flow, but not its fluctuations from event to event [14-16]. 
If one assumes that the elliptic flow is linearly propor- 
tional to the initial eccentricity, one can distribute the 
elliptic flow fluctuations around this average vi value in 
the same way as the initial eccentricity fluctuates around 
its average [17] and use this to predict flow fluctuations 
from eccentricity fluctuations [18]. This ignores, how- 
ever, the recently discovered [19] fact that initial-state 
shape fluctuations of the collision region lead to event- 
by-event fluctuations not only of the elliptic deformation 
e x , but simultaneously of all higher-order harmonic ec- 
centricity coefficients [20, 21], and that the simultane- 
ous presence of several harmonic eccentricity coefficients 
can lead to hydrodynamic cross-talk between anisotropic 
flows of different harmonic order [20] . Without the possi- 
bility to hydrodynamically evolve fluctuating initial con- 
ditions event-by-event, the assumption of a linear depen- 
dence of i>2 on e x can thus not be rigorously tested. 

Experimentally, anisotropic flow coefficients v n are 
measured by analyzing multiparticlc correlations in az- 
imuthal angle around the beam axis [24] . Evcnt-by-event 
flow fluctuations and additional non-flow correlations in- 
fluence the elliptic flow V2 derived from different such 
measures in different ways [22-25] and affect the magni- 
tude of the extracted flow. This has serious implications 
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for the determination of the quark-gluon plasma shear 
viscosity from elliptic flow data [5, 6]: quantitatively 
trustworthy results require a detailed understanding of 
the spectrum of fluctuations of the anisotropic flow coef- 
ficients in the experimentally observed final state and its 
relation to the fluctuations of the corresponding moments 
of the initial eccentricity distributions that are believed 
to drive the flows measured by various different meth- 
ods. In addition, widely held beliefs as to which mo- 
ments of the initial eccentricity distribution are directly 
related to which moments of the final flow distribution 
[17, 21, 26] must be tested by generating the final states 
hydrodyamically event by event, i.e. separately for each 
fluctuating initial condition [20, 27-31]. These are the 
goals addressed in the present article. (For related work 
using non-hydrodynamic models for the transfer of initial 
eccentricities into final flows sec [32].) 

To generate fluctuating initial conditions we use Monte 
Carlo versions of the Glauber [33] and f KLN [34, 35] mod- 
els, in the implementation by Hirano and Nara [11, 12]. 
After defining and compiling in Sec. II several differ- 
ent definitions of initial-state eccentricities and final-state 
harmonic flow coefficients, we explore in Sec. Ill the cen- 
trality dependence of the average ellipticity, triangular- 
ity and a few higher-order harmonic coefficients in their 
various incarnations found in the literature. In Sec. IV 
we use ideal event- by- event hydrodynamics 1 to analyze 
the correlations between final state flow anisotropies and 
their associated flow angles with the initial state ec- 
centricity coefficients and their associated angles. Wc 
identify strong cross-talk between coefficients of different 
harmonic order, especially in peripheral collisions with 
strong elliptic flow. In Section V we compare the conver- 
sion efficiency of initial-state ellipticity and triangularity 
into final-state elliptic and triangular flow in single-shot 
and event-by-event hydrodynamics. We find significant 
differences and discuss their implications. We summa- 
rize our results in Section VI and discuss different radial 
weights for the eccentricity definitions in the Appendix. 



II. DEFINITIONS 

In this section we discuss different definitions for the 
harmonic flow and eccentricity coefficients and briefly de- 
scribe the models used in computing the initial entropy 
and energy density profiles whose eccentricities are evalu- 
ated in Sec. Ill and which are evolved hydrodynamically 
in Sections IV and V. 



A. Ellipticity 

Usually known simply as "eccentricity" , we define the 
"ellipticity" e = £2 of a given matter distribution in the 
transverse (x, y) plane in terms of its r 2 -weighted second 
azimuthal moment [19, 26], 

£2 e l2 ^ F = I dxd y r2ei2<l ' e ( x ^y) 

/ dx dy r 2 e(x, y) 

where x = r cos <f>, y = r sin <f>. This formula assumes that 
the origin is the center of the distribution e(x,y). In a 
Monte Carlo approach for generating the initial distribu- 
tion e(x,y) (see Sec. II D) this must be ensured by re- 
centering each event before using Eq. (1). By default 
we characterize in Eq. (1) the matter distribution by 
its energy density e(x,y) [36]. Since some authors (e.g. 
[11, 12]) prefer defining the source ellipticity in terms of 
its entropy density distribution s(x, y), we compare in 
Sec. Ill B energy- and entropy- weighted ellipticitics. 

In Eq. (1), x and y are "reaction plane" (RP) coor- 
dinates: The reaction plane is the (x, z) plane, with z 
pointing along the beam and x pointing along the di- 
rection of the impact parameter b between the colliding 
nuclei, y is perpendicular to the reaction plane. Because 
of the minus sign on the r.h.s. of Eq. (1), the angle ip2 
on the l.h.s. of Eq. (1) points in the direction of the 
minor axis of the corresponding ellipse. For an cllipti- 
cally deformed Gaussian density distribution, this is the 
direction of the largest density gradient and thus of the 
largest hydrodynamic acceleration and also of the finally 
observed elliptic Sow. The direction of this minor axis 
defines, together with the beam direction z, the "partici- 
pant plane" (PP). It is tilted relative to the reaction plane 
by i{>2 P ■ The label "participant" is motivated by the fact 
that the initial energy and entropy density distributions 
of the collision fireball reflect (more or less directly, de- 
pending on the model for secondary particle creation) 
the transverse distribution of the nucleons participating 
in the particle production process. The ellipticity £2 in 
Eq. (1) is correspondingly called "participant eccentric- 
ity" and also denoted as £ pa rt- 2 It can be written as 



yj{r 2 cos^)} + {r 2 sin(2(^)} 
FT 

y/{y 2 -x 2 } 2 +4{xy} 2 
{y 2 +x 2 } 



For convenience we here use the longitudinally boost-invariant 
(2+l)-dimcnsional viscous hydrodynamic code VISH2+1 [8] for 
zero viscosity, restricting our analysis to anisotropic transverse 
flow at midrapidity. This allows for easy later inclusion of vis- 
cous effects. A (3+l)-d viscous hydrodynamic code has recently 
become available [30], and a comparison between (2+l)-d and 
(3+l)-d viscous evolution near midrapidity is in progress [48]. 



2 Traditionally e pa rt is defined in terms of the transverse density 
of wounded nucleons, but since what matters for the subsequent 
hydrodynamic evolution is not the distribution of wounded nu- 
cleons themselves but of the matter generated by the wounded 
nucleons, we use the name e pa rt for the ellipticity characterizing 
the thcrmalizcd matter. 
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Here {...}= J dxdy (. . . ) e(x, y) defines the "event aver- 
age" over the matter distribution e(x, y) in a single colli- 
sion event [23] . Equivalently, the participant eccentricity 
can be written as 



{y 2 -x 2 Y 
{y 2 +x 2 }' 



-part 



(3) 



where {...}'= J dxdy (...) e'(x,y) indicates the aver- 
age over a rotated event with energy density e'(x,y) = 
e(xcos'02 p — ysin , !/'2 P ,£sin'!/'2 P +y C0S V'2' P ) whose mi- 
nor and major axes now align with x and y. 

The event-average {. . . } is to be distinguished from the 

the "ensemble average" {•••) = jj^2n = ii- ■ • }« wnere N 
is the total number of events and {...}„ is the event- 
average over the energy density e n (x, y) in event number 
n. The average participant eccentricity is thus defined as 



1 



/v 



-part) — / J (£part)n- 

n = 1 



(4) 



This differs from the mean eccentricity £ par t of the av- 
erage (recentered and rotated by ip2 P ) energy density 
e'(x, y) = jj J2 n = i e n( x > v) which can be written in the 
following equivalent ways: 



-part 



y/{{y 2 -x 2 }) 2 +*{{xy}) 2 _ <{y 2 -* 2 }') 



{{y 2 +x 2 }) 



{{y 2 +x 2 }') ■ 



(5) 



by 



£rp 



{y 2 -* 2 } 

{y 2 +x 2 } 



(6) 



in terms of an event-average over the (properly cen- 
tered) energy density e(x,y). The so-called standard 
eccentricity is defined as the analogous ratio of expec- 
tation values taken with a smooth average energy den- 
sity e(x, y) = -h Sn = l e «( x > v) obtained by superimpos- 
ing many events without rotating them from the partici- 
pant to the reaction plane: 



£rp 



({y 2 -* 2 }) 

({y 2 +x 2 })' 



(7) 



In other words, the standard eccentricity is the mean 
reaction plane eccentricity. In contrast, the average re- 
action plane eccentricity is defined by 



(erp) = 



{y 2 



-x 2 } 



{y 2 +x 2 } 



(8) 



In contrast to (4), one here ensemble- averages over nu- 
merator and denominator separately before forming the 
ratio. 

Nature performs heavy-ion collisions event by event, 
and hydrodynamic forces generate in each event an ellip- 
tic component Vi of the anisotropic flow which is causally 
related to the specific initial ellipticity e pa rt in that event. 
Theorists often do not compute the hydrodynamic evolu- 
tion of the collision fireball event by event, but approxi- 
mate Nature's procedure by generating from a superposi- 
tion of many fluctuating initial conditions a single smooth 
initial distribution e(x 7 y) which they then evolve hydro- 
dynamically in a "single shot" , extracting the mean ellip- 
tic flow V2 corresponding to the mean eccentricity £ par t of an d 
that averaged source distribution. Obviously, V2 is a de- 
terministic consequence of e par t and does not fluctuate at 
all; it can not be measured experimentally. What can (at 
least in principle, although not easily in practice) be mea- 
sured experimentally [16] is the average elliptic flow (i^) 
of a large ensemble of collision events. This observable is 
conceptually more closely related to (e pa rt) than to £ pa rt; 
for an exactly linear hydrodynamic response V2 ~ £ pa rt, 
one has (v2)/(e P art) = ^/^part [17]. We will explore the 
differences between £ par t and (e par t) and discuss conse- 
quences for the theoretically computed V2 as opposed to 
the measured [16] (1)2) in Sees. Ill A and V. 

In addition to these "participant eccentricities" one 
can also define "reaction plane eccentricities" . For a sin- 
gle event, the reaction plane eccentricity £r P is defined 



Contrary to what the reader may have been led to be- 
lieve by our remarks above, experiments do not directly 
measure the average elliptic flow (772) (which for linear 
t>2 ~ £part would be directly related to the average par- 
ticipant eccentricity (£ par t) (3))- Instead they measure 
quantities such as i>2{EP}, i>2{2}, and ^{4} that, even if 
so-called non-flow contributions could be completely ig- 
nored, are affected by event-by-event ^-fluctuations and 
thus differ from (1)2) • (^2) can be reconstructed from 
the experimental measurements with some additional as- 
sumptions [16] which on the surface look harmless but 
should be further tested. Motivated by the hypothesis of 
linear hydrodynamic response, t>2 ~ £ pa rt, these V2 mea- 
sures motivate the definition of corresponding ellipticity 
measures [17], the so-called 2 nd and 4 th order cumulants: 



^{2} = \(e 2 t ) 



^{4}=((4art) 2 -((4art)-(4art) 2 )) 1/4 



(9) 



(10) 



Note that the last expression involves the difference of 
two positive definite quantities which itself does not need 
to be positive definite. If fluctuations get large, the ex- 
pression under the fourth root can become negative, leav- 
ing £{4} undefined. We will sec that this can happen in 
the most central and the most peripheral ccntrality bins. 

It was shown in [22] that in the MC-Glauber model the 
real and imaginary parts of the complex ellipticity defined 
by Eq. (1), with the wounded nucleon density as weight 
function on the r.h.s., both have approximately Gaussian 
fluctuations, with equal widths a e . If this is the case, the 
magnitude £2 of this ellipticity exhibits fluctuations of 
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Besscl-Gaussian type 3 [37], leading to the identity [22] 



; {4} 



£RP 



(11) 



For sufficiently large average ellipticities (£2) (i.e. suf- 
ficiently large impact parameters) one may hope to be 
able to ignore the restriction that e 2 can never fluctuate 
to negative values, and correspondingly assume the £2 
exhibits Gaussian (instead of Bessel-Gaussian) fluctua- 
tions. In this case one has [22] 



£{2} 2 = <e part ) 
£{4} 2 = 



2 1 al 



«£ par t} 2 - a 2 ) 2 - 2a 4 , 



(12) 



from which it follows that (e^ =£ P art) 4 is the arithmetic 
mean of e{2} 4 and e{4} 4 : 



£{2} 4 +£{4} 4 
2(£part) 4 



(13) 



We will use Eqs. (11) and (13) (which hold irrespective 
of the fluctuation width <r e ) in Sec. Ill A, and their ana- 
logues for the elliptic flow u 2 in Sec. V C, to test the as- 
sumptions of Bessel-Gaussian and Gaussian fluctuations 
of the event-by-event ellipticity and elliptic flow fluctu- 
ations in the Monte Carlo Glauber (MC-Glaubcr) and 
Monte Carlo fKLN (MC-KLN) models. 

If the hydrodynamic response were indeed linear, V2 ~ 
e par t , and non-flow effects could be ignored, the following 
identities would hold: 



(V2) V2 _ ^{2} _ « 2 {4} 

(^part) £part 

~ £{2} £{4} ■ 



(14) 



To test these theoretically one needs event-by-event hy- 
drodynamics which is the only possibility to properly ac- 
count for event-by-event flow fluctuations. In the past, 
event-by-event hydrodynamical evolution of fluctuating 
initial conditions has been technologically out of reach. 
Comparisons between theory and experiment have been 
based on "single-shot hydrodynamic evolution" which 
propagates a smooth initial profile obtained by cither 
using the so-called optical versions of the Glauber and 
fKLN models or averaging over many fluctuating initial 
profiles from their Monte Carlo versions (MC-Glauber 
and MC-KLN, respectively). Assuming linear hydrody- 
namic response, one can still compare the theoretically 
computed (1)2) with the experimentally measured i>2{2} 
or u 2 {4} if one normalizes the former by (e par t) an d the 
latter by £{2} or £{4}, respectively, calculated from the 
same initial state model [5, 6]. In this context the iden- 
tity £{4} = (erp) (which holds if the ellipticity fluctua- 
tions are Gaussian) becomes particularly useful because 



3 This takes into account that £2 can never fluctuate to negative 
values. 



it suggests that the measured U2{4} can be directly com- 
pared with a single-shot hydrodynamic V2 obtained from 
a smooth reaction-plane averaged initial density of ellip- 
ticity (£rp), without any corrections for flow fluctuations. 
Even better, 1*2 {4} can be shown to be completely free of 
two-particle non-flow contributions [17, 22]. These argu- 
ments have been used in [12] and provide a strong motiva- 
tion for us to test the underlying assumptions (Gaussian 
ellipticity fluctuations and linear hydrodynamic elliptic 
flow response) in the present work. 

We close this subsection by recalling the expression for 
the participant plane angle of a given event (see e.g. [23]) 



pp 
2 



■ tan 



2{xy} 
{y 2 — x 2 } 



and for its transverse area 



s = TrV{x 2 Y{y 2 Y- 



(15) 



(16) 



Both expressions assume that the events are properly 
centered at the origin. 



B. Higher order eccentricity coefficients 

The definition (1) can be generalized to higher har- 
monic eccentricity coefficients [19, 26]: 



£ n e 



inip T 



J dxdyr 2 e m ^ e(x,y) 
J dx dy r 2 e(x, y) 



(17) 



Alternatively one can use r n instead of r 2 as radial weight 
on the right hand side [20]: 



/ dxdy r n e mrt> e{x,y) 
J dx dy r n e(x, y) 



(18) 



Still another variant uses the entropy density s{x, y) in- 
stead of the energy density e(x, y) as weight function: 



£n(s)e m < P(s) = 
<(*)e* 



J dx dy r 2 e in ^ s(x,y) 
J dx dy r 2 s(x, y) ' 

/ dx dy rV^ s(x,y) 
J dx dy r n s(x, y) 



(19) 
(20) 



We note that the r 2 -weighted eccentricity coefficients e n 
fall off faster with increasing harmonic order n than the 
r n -weighted eccentricities e' n (see Appendix). Also, as in 
(1), the minus sign in Eqs. (17)-(20) guarantees that, for 
a Gaussian density distribution that has only n th -order 
eccentricity e n , the angle ■0^ p points in the direction of 
the steepest density gradient, and thus in the direction of 
the corresponding hydrodynamically generated n th -order 
harmonic flow v n (see next subsection). It can be written 
as — l = e _m W n ) and amounts to a rotation of tp^ by 
ir/n. For example, if the profile is square-shaped, ip, 
points to the sides instead of its corners. 



pp 
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A complete characterization of the fluctuating initial 
density profile, that captures all aspects of the location 
of "hot-spots" and their gradients, uses an expansion of 
the initial (energy or entropy) density profile in terms of 
cumulants [38]. We will postpone their discussion to a 
future analysis. 

As stated before, we will use the energy density as the 
default weight function; in cases of possible ambiguity, 
we will use the notations e„(e), e n (s) etc. to distinguish 
between energy and entropy density weighted eccentricity 
coefficients and angles. Eccentricities e without harmonic 
index n denote ellipicities (i.e. in the absence of n, n = 2 
is implied). 

The coefficients e n and angles ip pp define the eccen- 
tricies and angles of the matter distribution in the par- 
ticipant plane. We note that the participant plane angles 
tp pp associated with eccentricity coefficients of different 
harmonic order n do not, in general, agree (see Sec. IV A). 
We will not study higher harmonic generalizations of the 
reaction-plane ellipticity (6). 



C. Harmonic flow coefficients 

We charaterize the finally observed momentum dis- 
tribution dN / {dy prdpT d(j)p) by "harmonic flow coeffi- 
cients" constructed in analogy to Eq. (17), but without 
the extra minus sign: 



m0 p 

i o 



dN 



v n (y,p T )e in ^^ = — JvPTdvrd*^ (21) 



dyp T dp T 



v n (y) e 



dN 
dy 



(22) 



In boost-invariant hydrodynamics they are rapidity-inde- 
pendent, so we drop the argument y and keep in mind 
that we should only compare with midrapidity data at 
y = where the assumption of boost-invariant longitudi- 
nal expansion is most justified. The spectra -3 d , N 

1 ■> r dyp T dp T d<p v 

are computed from the hydrodynamic output with the 
Cooper-Frye prescription [39] along an isothermal ki- 
netic decoupling surface of temperature Tdcc = 140 MeV. 
Equation (21) defines the pT-differential harmonic flow 
Vu{pt) and flow angle 4>n P (PT), whereas Eq. (22) gives 
their ^-integrated values v n and tp^. The orientation 
of the final momentum distribution defines the "event 
plane" , indicated by superscript EP. Again, different har- 
monic flows are usually associated with differently ori- 
ented event planes. The first three harmonic flow coef- 
ficients are the directed flow (vi), elliptic flow (^2), and 
triangular flow (113). 



D. Initial-state models 

We use Monte Carlo versions [11, 12] of the Glauber 
[33] and f KLN [35] models to generate fluctuating initial 



conditions for the entropy density in 200 A GeV Au+Au 
collisions. For the MC-Glauber model we assume a two- 
component (soft+hard) model with a small hard frac- 
tion ((5 = 0.14 [11]); we also use a Woods-Saxon profile 
for the distribution of nucleon centers whose radius and 
surface thickness parameters have been corrected for the 
finite nucleon size [11]. The resulting entropy density 
profile is normalized to the final charged hadron multi- 
plicity density dN c \Jdy in central collisions; after this 
normalization, the centrality dependence of the initial 
entropy production is fixed by the model (MC-Glauber 
or MC-KLN). To convert the initial entropy density to 
energy density, we use the equation of state (EOS) s95p- 
PCE which matches Lattice QCD data at high tempera- 
tures to a chemically frozen hadron resonance gas at low 
temperatures [40, 41], using T c hom = 165 MeV as chemical 
freeze-out temperature. 



TABLE I: Centrality table for Au+Au at 200,4 GeV [11]. 



centrality 


&min(frn) 


bmax(fin) 


b (fm) 


Apart 


0-5% 


0.0 


3.3 


2.2 


352.2 


5-10% 


3.3 


4.7 


4.04 


294.7 


10-15% 


4.7 


5.8 


5.27 


245.6 


15-20% 


5.8 


6.7 


6.26 


204.2 


20-30% 


6.7 


8.2 


7.48 


154.5 


30-40% 


8.2 


9.4 


8.81 


103.8 


40-50% 


9.4 


10.6 


10.01 


64.9 


50-60% 


10.6 


11.6 


11.11 


36.6 


60-70% 


11.6 


12.5 


12.06 


18.8 


70-80% 


12.5 


13.4 


12.96 


7.5 


80-90% 


13.4 


14.3 


13.85 


4.4 



In the following we compute harmonic eccentricity and 
flow coefficients as functions of impact parameter b and 
collision centrality (%). The centrality classes are de- 
fined in terms of percentages of the total inelastic cross 
section, calculated from the distribution of the num- 
ber of wounded nucleons dN eveak / dN pait in the optical 
Glauber model (i.e. without accounting for fluctuations 
in Vp ar t at given impact parameter). Each centrality 
class is thus characterized by a range of impact param- 
eters & m i n < b < & max and an average value b, together 
with a mean number of wounded nucleons A par t. They 
are listed in Table I [11]. 



E. Averaging procedures for the initial profiles 

In this work we compare results obtained from an 
event- by- event hydrodynamical evolution of fluctuating 
initial conditions with the traditional method of "single- 
shot" hydrodynamic evolution, where one first averages 
over many fluctuating initial profiles to obtain a smooth 
average profile, and then evolves this smooth profile hy- 
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20 40 60 80 20 40 60 80 

Centrality (%) Centrality (%) 



FIG. 1: (Color online) Different ellipticities as a function of impact parameter (top row) or collision centrality (bottom row), 
for the MC-Glauber (panels (a) and (c)) and the MC-KLN model (panels (b) and (d)). Panels (a,b) show e-weighted, panels 
(c,d) show s-weighted ellipticities. (See Figs. 3, 4 below for a direct comparison between e- and s-weighted eccentricities.) Open 
stars indicate negative values for e{4} 4 . 



drodynamically. The question addressed in this compar- 
ison is to what extent the average harmonic flow coeffi- 
cients from event-by-event hydrodynamics can be faith- 
fully represented by the harmonic flow coefficients ex- 
tracted (at much lower numerical expense) from the hy- 
drodynamic evolution of an "average event" . 

Taking the initial density profiles from the Monte Carlo 
generator and superimposing them directly without ad- 
ditional manipulations generates a "reaction plane av- 
eraged" profile with ellipticity £rp (Eq. (7)). After re- 
centering each event to the origin of the x-y-planc, we 
can compute event by event the reaction and participant 
plane ellipticities (Eqs. (6) and (2,3)) and evaluate their 
ensemble averages (8) and (4), respectively. To gener- 
ate a smooth average profile with ellipticity e pa rt (Eq. 5) 
we rotate each recentered event by the angle ip2 F ( e ) 
(ip^i 3 )) ^ we want to determine the eccentricity of the 
average energy (entropy) density. For the calculation of 
entropy-weighted average eccentricities we perform any 
ensemble average first and convert the result to energy 



density later; in this case all events are rotated by their 
-02 P ( S ) angles. For energy- weighted ensemble averages 
we convert s to e in each event first, rotate by ip2 P ( e ) 
and perform the ensemble average last. Other sequences 
or mixtures of these steps are technically possible but 
physically not meaningful. Note that the processes of 
computing the energy density from the entropy density 
via the EOS and of averaging the event profiles do not 
commute: The energy density obtained via the EOS from 
the ensemble-averaged entropy density profile is not the 
same as the ensemble-averaged energy density where the 
EOS is used in each event to convert s to e. 

III. ECCENTRICITIES 

A. Centrality dependence of different ellipticities 

Fig. 1 shows a comparison between the different ellip- 
ticities defined in Sec. II A, as functions of the impact 
parameter b in panels (a) and (b) and as functions of col- 
lision centrality (as denned in Sec II D) in panels (c) and 
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(d). For panels (a) and (b) we generated 10,000 initial 
profiles for each impact parameter (except for 6 = 0, 1, 
and 2fm for which we generated 30,000 events each); for 
panels (c) and (d) we averaged over 10,000 profiles for 
each centrality bin. Within the centrality bins, the im- 
pact parameters were sampled between 6 m i n and & max 
with bdb weight. Compared to panels (a) and (b), this 
leads to additional cllipticity fluctuations related to the 
fluctuating impact parameter, whereas in Fig. la,b only 
-Wpart fluctuations at fixed b contribute. 

As discussed in Sec. IIA, Eq. (10), e{4} 4 can become 
negative when fluctuations grow large. Whenever this 
happens, we replace e{4} by \/|e{4} 4 | and indicate this 
by an open star in Fig. 1 (connected by dotted lines to 
other points in the graph) . One sees that e{4} 4 has a ten- 
dency to turn negative in the most peripheral collisions. 
In very central collisions ^{4} becomes very small, with 
central values that can have either sign depending on 
whether we keep the impact parameter fixed (Figs. la,b) 
or average over events with different impact parameters 
in a given centrality bin (see the 0—5% centrality values 
in Figs. lc,d). Statistical errors are large, however, and 
within errors e{4} 4 is compatible with zero for impact 
parameters b < 3fm, i.e. in the most central (0—5% cen- 
trality) collisions. We indicate this by open-ended error 
bars for -^/|e{4} 4 | , pointing from its upper limit all the 
way to zero. 

Comparing panels (a,c) for the MC-Glaubcr model 
with panels (b,d) for the MC-KLN model we see great 
similarities in shape, but systematic differences in mag- 
nitude of the ellipticities. The ratio of the MC-KLN and 
MC-Glaubcr ellipticities is shown in Fig. 2. Except for 
the most central and most peripheral collisions, the MC- 
KLN ellipticities exceed the MC-Glauber ones by an ap- 
proximately constant factor of around 1.2. Please note 




6 8 

b (fin) 



FIG. 2: (Color online) Impact parameter dependence of the 
ratio of ellipticities obtained from the MC-KLN and MC- 
Glauber models as shown in Figs. la,b. 



the difference in the ratios for the reaction plane and 
participant eccentricities at small b. (The point for £rp 
at b = 1 fm is obtained from a ratio of very small num- 
bers and probably not statistically robust - we had only 
30,000 events to determine the ensemble-averaged den- 
sity profile.) For the e{4} ratio we dropped all points 
where the statistical error for e{4} 4 extended into the 
region of negative values for either the MC-Glauber or 
MC-KLN model. 

Figure 1 shows that, for central and mid-peripheral 
collisions, the ensemble-averaged participant and reac- 
tion plane eccentricities (e pa rt) and (erp) agree very well 
with the mean eccentricities e par t and epp of the cor- 
responding ensemble-averaged profiles. For strongly pe- 
ripheral collisions (6 > 10 fm), however, the average of 
the ratio (Eqs. (3,4,8)) differs strongly from the ratio 
of averages (Eqs. (5,7)), indicating strong cvent-by-event 
fluctuations. We note that in very peripheral collisions 
the average event ellipticity drops quickly with increasing 
b while the cllipticity of the average profile remains large; 
single-shot hydrodynamic calculations based on a smooth 
average initial profile thus overestimate the effective ini- 
tial source ellipticity and produce more elliptic flow than 
expected from event-by-event hydrodynamic evolution of 
individual peripheral events. Still, as first emphasized in 
[6] , the calculated v-i from single-shot hydrodynamics de- 
creases steeply at large collision centralities [12, 42-44], 
due to the decreasing fireball lifetime, which contrasts 
with the initially reported experimentally observed be- 
haviour [45, 46], but agrees qualitatively with a recent 
reanalysis [47] where non-flow effects have been largely 
eliminated and/or corrected for. We do point out that 
our Monte-Carlo simulations do not include fluctuations 
in the amount of entropy generated per nucleon-nucleon 
collision [20]; these could have important effects on the 
ellipticities in very peripheral collisions. 

Comparing the curves for (e pa rt)j ^{2} and e{4} in 
Fig. 1 we see that (as is manifest in the Gaussian model 
analysis in Eq. (12)) e{2} receives a positive and £{4} 
receives a negative contribution from event- by- event el- 
lipticity fluctuations. In Fig. 3 we check, as a func- 
tion of impact parameter, the validity of the identities 
(11) and (13) which follow from Bessel-Gaussian and 
Gaussian e par t distributions, respectively. We see that 
both hold with good accuracy in the mid-centrality range 
(6<10fm for Eq. (11), 5<6<llfm for Eq. (13)) but 
break down in the most peripheral collisions. Both the 
Gaussian and Bessel-Gaussian hypotheses work slightly 
better for the MC-KLN than for the MC-Glauber model. 
Consistent with the analysis in Ref. [22], the Gaussian 
fluctuation hypothesis for e par t breaks down at small 
impact parameters whereas (as theoretically expected 
[22]) the Bessel-Gaussian hypothesis appears to con- 
tinue to hold, although we are unable to make this 
statement with statistical confidence. (For the ratio 
(e{4} — (eRp)/(e{4}+(£'Rp) we again dropped all points 
for which the error band for e{4} 4 reaches into negative 
territory.) 
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FIG. 3: (Color online) Checks of the assumption of Gaus- 
sian and Bessel-Gaussian fluctuations for e pa rt (see text for 
discussion). 




FIG. 4: (Color online) e par t(e) vs. e pa rt(s) for 6000 initial 
profiles from the MC-KLN model (1000 each for & = 0, 0-5%, 
15-20%, 20-30%, 30-40%, and 50-60% centrality). 



The assumption of Gaussian fluctuations of the real 
and imaginary parts of the complex ellipticity (1) is of- 
ten used to argue that the average reaction-plane ellip- 
ticity (erp) can serve as a proxy for e{4} (see Eq. (11)), 
and that therefore reaction-plane averaged initial den- 
sity profiles can be used in single-shot hydrodynamics 
(which ignores event-by-event fluctuations) to simulate 
the experimentally measured t>2{4} values. Fig. 1 and 
the bottom curves in Fig. 3 show that i>2{4} values ob- 
tained from single-shot hydrodynamic simulations with 
reaction-plane averaged initial conditions [12, 43] should 
not be trusted quantitatively for centralities > 40%. 

To summarize this subsection, all the simplifying as- 
sumptions that allow to focus attention on the three 
quantities (£ par t), e{2} and e{4} only (by substituting 
(e P art) for £part and e{4} for e RP or (e RP )) hold well for 
central to mid-central collisions ( < 40% centrality) but 
break down for peripheral collisions. For > 40% central- 
ity there exists no substitute for event-by-event hydro- 
dynamics if one aims for quantitative precision in the 
comparison with experimental data, since the latter are 
strongly affected by non-Gaussian event-by-event fluctu- 
ations at those centralities. 



B. Ellipticities with different weight functions 

Figure 4 shows a comparison between the energy- and 
entropy- weighted ellipticities of the initial profiles gener- 
ated with the MC-KLN model, on an event-by-event ba- 
sis. The scatter plot is based on 6000 events, 1000 each 
for b = and for the following finite-width centrality 
bins: 0-5%, 15-20%, 20-30%, 30-40%, and 50-60%. 
This is not a realistic mix in the experimental sense, but 
permits us to explore the full range from very small to 
very large event ellipticities. The blue dots in Fig. 4 rep- 



resent bin averages, and the solid black line is a linear fit 
through the origin. The fitted slope is 1.00, the scatter 
plot is seen to be tightly clustered around this fitted line, 
and only at small ellipticities £2 < 20% the e- weighted 
values are seen to be slightly larger on average than their 
s- weighted counterparts (see also Fig. 5a below). 



C. Higher order harmonics 

In Figs. 5a-d we compare the centrality dependences 
of the ensemble-averaged second to fifth harmonic eccen- 
tricity coefficients (energy- and entropy-weighted) from 
the MC-Glauber and MC-KLN models. The contour 
plots give a visual impression of the degree of deformation 
corresponding to the (larger) MC-KLN eccentricities, as- 
suming (for illustration) the absence of any other eccen- 
tricity coefficients than the one shown in the particular 
panel. 

First, one observes very little difference between the 
eccentricities of the entropy and energy density profiles, 
except for very central collisions (b < 5 fm for the MC- 
Glauber, 6<3fm for the MC-KLN model) where the 
energy-weighted eccentricities lie systematically some- 
what above the entropy-weighted ones (for all orders n 
studied here). The difference between s- and e-weightcd 
eccentricities at small b is bigger in the MC-Glauber than 
in the MC-KLN model. 

Next, one notes the significantly larger ellipticities 
and quadrangularities of the MC-KLN distributions com- 
pared to those from the MC-Glauber model, for all but 
the most central collisions. These are driven by geome- 
try, i.e. by the almond-shaped deformation of the nuclear 
overlap zone in non-central collisions, which in the KLN- 
modcl is more eccentric than in the Glauber model. The 
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FIG. 5: (Color online) Harmonic eccentricity coefficients £2 =£ P art (a), £3 (b), £4 (c) and £4 (d) as functions of impact parameter, 
calculated from the MC-Glauber (filled symbols, solid lines) and MC-KLN models (open symbols, dashed lines), using the 
energy density (circles) or entropy density (triangles) as weight function. The contour plots illustrate deformed Gaussian 

profiles e(r,(f>) = eg exp j^— (l+e n cos(n^>))J , with eccentricity e n (e) taken from the MC-KLN model at the corrsponding 
impact parameter. 




FIG. 6: (Color online) Scatter plots illustrating the event-by-event correlation of £3 (a), £4 (b), and £5 (c) with the fireball 
ellipticity £2 = £ par t , for the same set of 6000 event profiles analyzed in Fig. 4. The blue dots connected by a line indicate £2-bin 
averages, to guide the eye. 
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third and fifth order harmonics, which arc entirely due to 
fluctuations (and whose associated angles ip pp are there- 
fore completely uncorrclated to the reaction plane - see 
Rcf. [20] and discussion below), show remarkably simi- 
lar eccentricity values in the two initialization models, 
except for the most peripheral events. Comparing the 
viscous suppression of elliptic and triangular flow thus 
should allow to distinguish experimentally between the 
MC-Glauber and MC-KLN models [49]. 

Third, in central collisions all four eccentricity coeffi- 
cients are roughly of the same size. In peripheral colli- 
sions, the fluctuation-dominated eccentricity coefficients 
(e 3 and e 5 ) are generically smaller than the geometry- 
dominated ones (e 2 , but also to some extent e 4 ). 4 This 
is less obvious when one defines the higher order eccen- 
tricities with r" instead of r 2 weight [20] , which tends to 
increase the values of the higher harmonics in peripheral 
collisions. 

Even with "only" an r 2 weight, £4 and £5 are seen 
to become large enough around b ~ 10 — 13fm that, 
if collective acceleration happens predominantly in the 
directions of steepest descent of the density profile, one 
has to expect cross-currents in the developing anisotropic 
flow patterns. These can lead to destructive interference 
and a correspondingly reduced efficiency of converting 
n th -order eccentricities e n into n th -order harmonic flows 
v n [26]. In realistic situations this issue is exacerbated 
by the simultaneous presence of several large eccentricity 
components e n , which is expected to lead to a strongly 
non-diagonal and probably non-linear response matrix 
relating v n to e n [20]. This will be discussed in Sec. IV. 



D. Eccentricity correlations 

It is reasonable to ask whether and how the different 
harmonic eccentricity coefficients e n are correlated with 
each other. Figure 6 shows scatter plots of the correla- 
tions between £3.4,5 and the cllipticity £2 which, for large 
£2 values, is dominated by geometric overlap effects. We 
note that, according to the definition (17), all eccentric- 
ity coefficients are positive definite, e n > 0. Keeping this 
in mind, Figs. 6a, c show that £3 and £5 are uncorrelated 
with the fireball ellipticity; the slight growth of (£3,5) 
with increasing £2 is related to the growth of the vari- 
ances of their distributions in more peripheral collisions. 

In contrast, the quadrangularity £4 shows a clear pos- 
itive correlation with the ellipticity, see Fig. 6b. It is 
of geometrical origin: it reflects the football or almond 
shape of the overlap zone in non-central collisions which 



4 We checked that the centrality dependences of the ratios e n /sz 
agree qualitatively, but not quantitatively with Fig. 3 in Rcf. [21]. 
We suspect that the differences, which are larger for the MC- 
Glauber than the MC-KLN model, are due to somewhat differ- 
ent Woods-Saxon and (in the MC-Glauber case) fluctuation size 
parameters used in Rcf. [21]. 



is a little sharper than a pure cos(2</>) deformation. This 
is corroborated by the angle ip P p shown in Fig. 7a which, 
on average, points 45° relative to ip PP (which again points 
in x-direction) . This means that the quadrangular com- 
ponent of the initial fireball definition is oriented like a 
diamond, with its corners on the x and y axes. Super- 
imposing it on a pure cos(2(/>) deformation leads to a 
somewhat sharper shape of the density distribution. 



IV. EVENT-BY-EVENT HYDRODYNAMICS 
AND FLOW FLUCTUATIONS 

In this section we analyze the results from ideal fluid 
event-by-event hydrodynamic evolution of the fluctuat- 
ing initial profiles studied in the previous section. We 
focus on the anisotropic flow coefficients v n , their rela- 
tionship to the initial eccentricity coefficients e n , and the 
correlation between the n th -order flow angles i/j® p and 
the corresponding n th -order participant-plane angles tp pp 
associated with s n . 



A. Correlations between participant plane, event 
plane, and reaction plane 

One of the key characteristics of fluid dynamics is its 
ability to transform initial geometric deformation into a 
deformation of the final momentum distribution, via col- 
lective flow. This happens through spatially anisotropic 
hydrodynamic forces (i.e. pressure gradients) which 
cause anisotropic acceleration of the fluid. As a result, 
correlations between participant and event planes are ex- 
pected: The angle tpn F points in the direction of the 
largest pressure gradient associated with the n th har- 
monic component of the spatial deformation of the initial 
density distribution, while ip^ p points into the direction 
where the n th harmonic component of the final collective 
flow is largest. Without interference between harmonics 
of different order, we would thus expect , PP and tp^ p to 
point, on average and up to event-by-event fluctuations, 
in the same direction. 

In Figs. 7a, b we show the distribution of participant 
and event plane angles, associated with the rt th -order ec- 
centricities and harmonic flows, relative to the x-z reac- 
tion plane. The analysis uses the same 6000 events as 
before. In panel (a) we see that ^35 are completely un- 
corrclated with the reaction plane [20] , as expected from 
the fact that the corresponding eccentricities are entirely 
fluctuation-driven, without contribution from the colli- 
sion geometry Panel (b) shows that the same holds true 
for ipf p , which is (at least superficially) consistent with 
the expectation that ^3 is mostly or entirely driven by 
£3, and V5 by £5. We will revisit this below. ip PP and 
■0f p are strongly correlated with the reaction plane, at 
least for this mixed-centrality set of events. This is ex- 
pected since, for non-central collisions, £2 is mostly con- 
trolled by the almond-shaped overlap geometry, and v 2 
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FIG. 7: (Color online) Event-by-event correlation of the participant plane (PP,(a)) and event plane (EP,(b)) angles with the 
reaction plane (RP), as well as the correlation between participant and event plane angles (c), for different harmonic eccentricity 
and flow coefficients. The same 6000 events as in Fig. 4 were analyzed. 



is mostly a collective flow response to this geometric de- 
formation; event-by-event fluctuations contribute to £2 
(and thus V2), but in general do not dominate them. 

The behavior of p in Fig. 7a is interesting because 
it is on average strongly "anti-correlated" with the re- 
action plane, in the sense that it points (on average) 
at 45° relative to the x-axis. The geometric reason for 
this has already been discussed above in subsection III D. 
On the other hand, Fig. 7b shows that the angle tfjf F 
points on average into the reaction plane. This correla- 
tion of t/jf p with the reaction plane is somewhat weaker 
than the anti-correlation of ?/> PP with that plane seen in 
panel (a). Still, it suggests that quadrangular flow V4 
does not, on average, develop predominantly in the di- 
rection of the steepest pressure gradient associated with 
£4, but in the direction of steepest £2-induced pressure 
gradient. This can be understood as follows: since £2 
generates a second harmonic deformation of the flow ve- 
locity profile which elliptically deforms the exponent of 
the flow-boosted Boltzmann factor exp[— p ■ u(x)/T(x)] 
describing the local thermal momentum distribution of 
particles, it leads to harmonic contributions V2k of all 
even orders n = 2k in the momentum distributions of the 
finally emitted particles [50]. Fig. 7b suggests that, on 
average, this effect wins over initial-state quadrangular 
deformation effects. 

Figure 7c, however, in which we analyze directly the 
correlation between the event and participant plane an- 
gles, paints a more subtle picture. It shows, surprisingly, 
a correlation peak at zero relative angle between iftf p and 
i/> PP , whereas the above discussion should have led us to 
expect a correlation peak at 45°. The resolution of this 
paradox is presented in the next subsection: The relative 
importance of geometric and fluctuation-induced contri- 
butions to e n , v n , and their associated angles changes 
with collision centrality, with geometry playing a rela- 
tively larger role in peripheral collisions. One should 
therefore look at the angle correlations as a function of 
collision centrality. One finds that the correlation func- 
tion peaks in Figs. 7a, b for the 4 th -order angles relative 



to the reaction plane are almost entirely due to geometric 
effects in peripheral collisions, while in central collisions 
both ■i/j PP and are fluctuation-dominated and thus 
essentially uncorrelatcd with the reaction plane. On the 
other hand, precisely because in central collisions geomet- 
ric effects such as geometrically driven elliptic flow do not 
dominate the hydrodynamic response to the fluctuation- 
driven higher-order eccentricities, ijjf p and ^ PP remain 
relatively strongly correlated in near-central collisions. 
This is the reason for the peak at 0° for n = 4 in Fig. 7c. 
(A hint of the "anti-correlation" at 45° is still visible in 
Fig. 7c, and it would be stronger if we had not (for unre- 
lated reasons) strongly oversampled central collisions in 
our mixed-centrality sample.) 



We close this discussion with the following additional 
observations about Fig. 7c: (i) The second-order partici- 
pant and event planes are much more strongly correlated 
with each other than either one of them is with the re- 
action plane. This shows that even in very central colli- 
sions, where the source ellipticity is mostly fluctuation- 
driven and its angle therefore only weakly correlated with 
the reaction plane, elliptic flow develops event- by- event 
in the direction of the short axis of the ellipsoid, (ii) 
Even though the angles associated with £3 and V3 are 
uncorrelated with the reaction plane (Figs. 7a, b), they 
are strongly correlated with each other. This indicates 
that i>3 is mostly driven by £3, especially in the more 
central collisions, with relatively little interference from 
other harmonics, (iii) The 5 th -order event and partici- 
pant plane angles show correlation peaks both at and 
7r/5. As we will see in the following subsection, the for- 
mer results from central and the latter from peripheral 
collisions. The peak at 7r/5 indicates significant cross- 
feeding between modes with n =2, 3, and 5. 
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FIG. 8: (color online) Event-by-event correlation between the absolute values of the event plane (EP) and participant plane 
(PP) angles for the harmonics of order 2 — 5 (panels (a)-(d)), for events in different centrality classes as indicated in the legend. 
(The same set of MC-KLN events as in Fig. 4 was used.) The insets show the centrality dependence of the widths of these 
correlations around zero. 



B. Centrality dependence of event and participant 
plane correlations 

Figure 8 looks at the correlation between the n th - 
order EP and PP angles at different collision central- 
ities. This generalizes a similar analysis for n = 2 in 
Ref. [29] to higher harmonics. Plotted are the distri- 
butions of the absolute value of the difference between 
the two angles in the main graph and the rms of this dis- 
tribution (i.e. the width around zero of the correlation) 
in the inset, as a function of collision centrality. Panel 
(a) shows that the second-order participant and event 
planes are strongly correlated at all collision centralitics. 
This demonstrates that elliptic flow is generated almost 
exclusively by the source ellipticity. The variance of the 
correlation is ~ 0.05 rad in the mid-central range (15-40% 
centrality) and increases in very central and very periph- 
eral collisions due to growing ellipticity fluctuations. 

A similar correlation exists for the 3 rd -order partici- 



pant and event planes, at all collision centralitics, but 
with a larger variance of order 0.2 — 0.3 rad (depending 
on centrality). The relatively strong correlation suggests 
that £3 is the dominant driver for 1*3 [19]. 

For the 4 th - and 5 th -order participant and event planes 
the situation is complicated, as seen in panels (c) and (d). 
The planes are correlated with each other (i.e. the dis- 
tributions peak at zero difference angle) in central colli- 
sions, become essentially uncorrelatcd in mid-central col- 
lisions and anti-correlated (i.e. peaked at a difference 
angle of %/n, n = 4, 5) in peripheral collisions. The anti- 
correlation in peripheral collisions indicates strong mode- 
mixing, driven by the large ellipticity £2 and strong el- 
liptic flow vi at large impact parameters which generates 
i>4 and V5 contributions by coupling to lower harmonics, 
as described in the previous subsection. For V4 in par- 
ticular, a strong cos(20) component in the collective flow 
velocity generates a U4 of the final momentum distribu- 
tion, without any need for nonzero £4. At large impact 
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FIG. 9: (Color online) (v n )(e„) for 71 = 2,3,4,5 (panels (a)-(d)). As in Fig. 4, each centrality class contains 1000 MC-KLN 
events, grouped in eccentricity bins of sufficient width to have reasonable statistics in each bin. 



parameters, ^-induced quadrupolar flow from the initial 
elliptic deformation of the overlap region thus dominates 
over any contribution from initial quadrangular deforma- 
tion. In near-central collisions, on the other hand, where 
all £„ stem mostly from shape fluctuations, 1^5 are dom- 
inantly driven by £4^. 

C. Harmonic flows and their corresponding initial 
eccentricities: nonlinear hydrodynamic response 

As discussed in the Introduction, it is often assumed 
that the harmonic flows v n respond linearly to the eccen- 
tricities e„, at least as long as the latter are small. This 
assumption receives support from hydrodynamic simula- 
tions [26] as long as one probes deformed initial profiles 
with only a single non-vanishing harmonic eccentricity 
coefficient. In Fig. 9 we investigate the validity of this 
assumption with fluctuating MC-KLN events which fea- 
ture nonzero e n values for all n. 

Figure 9a generally provides support for the assump- 
tion of a linear dependence of the elliptic flow v 2 on initial 
cllipticity £2, with two important caveats: 



(i) At small and large ellipticities, vi deviates upward 
from a best-fit line through the origin, indicating 
additional contributors to the elliptic flow. Indeed, 
for zero ellipticity £2 = we find a nonzero average 
(^2}. These are events with typically large nonzero 
values for eccentricities of higher harmonic order 
which generate elliptic flow through mode-mixing 
(e.g. between £3 and £5). We see that this happens 
at all centralities, even for 6 = 0, due to event-by- 
event fluctuations of the eccentricity coefficients. 

(ii) The slope of the curve (f2)(£2) decreases in very 
peripheral collisions, indicating destructive inter- 
ference via mode-mixing from other harmonics in 
the hydrodynamic evolution of the small and highly 
fluctuating fireballs created at large impact param- 
eters. 

The £3-dependence of triangular flow (v^), shown in 
Fig. 9b, shows a qualitatively similar story, but the devi- 
ations from linear response are stronger, with significant 
non-zero triangular flow in events with zero initial trian- 
gularity, especially for larger impact parameters. 
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FIG. 10: (Color online) Transverse momenrum spectra for directly emitted ir + (a) and protons (b) from event-by-event (solid 
lines) and single-shot hydrodynamics (dashed lines), for 200 A GeV Au+Au collisions at five selected impact parameters. 



For (1)4) and (vs), shown in Figs. 9c and 9d, mode- 
mixing effects are very strong, and a linear response of 
v n to e n (n = 4, 5) can no longer be claimed. This 
is quite different from the results in [26] where U4 was 
studied for a source that had only £4 deformation: in 
this case v^e^) was found to be approximately linear for 
small £4, with a downward bend at larger £4 values due to 
negative interference from cross-currents for sources with 
large quadrangularities. (This approximately linear de- 
pendence survived in the pT-integrated i>4 even though it 
was noticed in a related study [51] that, for mid-central 
collisions, the differential quadrangular flow v±{pt) ap- 
pears at high pt to be mostly determined by the elliptic 
deformation of the hydrodynamic flow profile generated 
by £2-) Our present study shows that it is unlikely that 
the anisotropic flow resulting from highly inhomogeneous 
initial profiles with nonzero eccentricity coefficients of all 
harmonic orders can be obtained by some sort of lin- 
ear superposition of flows generated from sources with 
only a single nonzero harmonic eccentricity coefficient, 
as suggested [38]. The hydrodynamic response {«„} to 
a set of initial eccentricity coefficients {£„} is not only 
non-diagonal, but also (via mode-mixing) non-linear, and 
there is no suitable single-shot substitute for event-by- 
event hydrodynamic evolution of fluctuating initial con- 
ditions. 

Wc note, however, that non-linear mode-mixing effects 
appear to be minimal for the elliptic and triangular flow 
(Figs. 9a, b). vi and W3 remain therefore the best candi- 
dates for an extraction of the fluid's viscosity, by study- 
ing (with quantitative precision) the fluid's efficiency in 
converting initial spatial deformations into final momen- 
tum anisotropics and anisotropic flows. Wc will further 



elaborate on this theme in the next section. 



V. SINGLE-SHOT VERSUS EVENT-BY-EVENT 
HYDRODYNAMICS 

We now discuss the effects of event-by-event initial- 
state fluctuations on the finally observed pion and proton 
p-r-spectra and anisotropic flow, comparing traditional 
single-shot hydrodynamic evolution of an appropriately 
constructed smooth average initial profile with event-by- 
event evolution of fluctuating initial conditions (with an 
ensemble average taken at the end). Since the calcu- 
lation of resonance decay feeddown corrections is com- 
putationally expensive but not expected to cause quali- 
tative changes, we here concentrate on directly emitted 
("thermal") pions and protons. For the graphs shown 
in this section, we generated for each impact parameter 
1000 fluctuating events and propagated them either indi- 
vidually ( "event- by-event hydrodynamics" ) or in a single 
hydrodynamic run after rotating and averaging their pro- 
files ( "single-shot hydrodynamics" ) down to a decoupling 
temperature of 120 MeV. 



A. Transverse momentum spectra 

In a very interesting recent paper [31] Chatterjee et 
al. showed that thermal photon spectra from exploding 
heavy-ion collision fireballs with fluctuating initial con- 
ditions which were hydrodynamically evolved event-by- 
event are significantly harder than those obtained from 
single-shot hydrodynamic evolution of the corresponding 
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FIG. 11: (Color online) Eccentricity-scaled elliptic flow v^jsi as function of impact parameter, for pions (panels (a,b)) and 
protons (panels (c,d)), from the ideal fluid dynamic evolution of initial MC-Glauber (a,c) and MC-KLN (b,d) density profiles. 
Solid (dashed) lines correspond to event-by-event (single-shot) hydrodynamics. See text for discussion. 



cnscmblc-avcragcd (and therefore much smoother) ini- 
tial profiles. The authors of [31] attributed this effect 
to the existence of "hot spots" in the fluctuating initial 
conditions that radiate photons at a higher than average 
temperature. Figure 10 shows that the same hardening 
occurs in the pion and proton spectra even though these 
strongly interacting hadrons are emitted only at freeze- 
out, with the same decoupling temperature assumed in 
both types of evolution. 5 This proves that the effect is 
due to stronger radial flow in the event- by- event evolved 
fluctuating fireballs, driven by the stronger than average 



5 A similar effect was also seen in [29] whose authors further 
pointed out that the strength of this "hardening effect" depends 
on the fluctuation size parameter in the initial conditions (i.e. 
the area over which the entropy produced in a nucleon-nucleon 
collision is distributed). 



pressure gradients associated with the "hot spots" (i.e. 
over-dense regions) in the initial profile. The importance 
of initial-state fluctuation effects on the final pp-spectra 
becomes stronger in peripheral collisions where the ini- 
tial fireballs are smaller and "hot spots" have a relatively 
larger influence. If stronger radial flow is the explana- 
tion of the fluctuation-driven hardening of the pion and 
proton spectra observed in Fig. 10, it is probably also 
a dominant contributor to the hardening of the photon 
spectra noted in Ref. [31], at least for low pr (i.e. in 
the hydrodynamic regime). This could be checked by 
comparing the photon radiation from the late hadronic 
stage in event-by-event and single-shot hydrodynamics 
which, if our interpretation is correct, should show the 
same fluctuation-driven, flow-induced hardening as the 
total photon spectra. 



16 



0.2 



0.15 



1U 



0.1 



0.05 





0.35 




0.3 




0.25 


o 
o 






0.2 


ro 




u 


0.15 


r- 






0.1 




0.05 



(a) 

5 *> MC-Glauber 




- + -v 3 /e 3 (e) (^ PP (e)-rot. & avg v 

- - -v 3 /e 3 (s) (^ p (e)-rot. & avg 
-*-v 3 /e 3 (s) (■0 PP (s)-rot. & avg.) 
+ (v 3 )/(e 3 (e)) 

- <i*)/<e 3 (<0> 
■ — iv3)/e 3 (e) 







6 8 

6 (fin) 



10 12 14 



MC-Glauber 




(c). 



6 8 10 12 14 

b (fin) 



0.2 



0.15 



10 

0.1 



0.05 


0.4r 

0.35 
0.3 
¥ 0.25 

O 

3 0.2 

CO 

10 

^0.15 
0.1 
0.05 

0- 





(b) 

MC-KLN 






7T+ 


V 


e-by-e hydro 

single-shot hydro 


V- - 



6 8 10 12 

b (fm) 



14 



(d) 



MC-KLN 




6 8 

b (fin) 



10 12 14 



FIG. 12: (Color online) Same as Fig. 11, but for the eccentricity-scaled triangular flow V3/S3. 



B. Elliptic and triangular flow 



In Figures 11 and 12, we compare the eccentricity- 
scaled elliptic and triangular flows, v 2 /s2 and v 3 /e 3 , for 
pions and protons as a function of impact parameter, 
from single-shot (dashed lines) and event-by-event hy- 
drodynamics (solid lines). These ratios represent the ef- 
ficiency of the fluid for converting initial spatial defor- 
mations into final-state momentum anisotropics. This 
conversion efficiency is affected (i.e. reduced) by shear 
viscosity, so these ratios form the basis of many anal- 
yses that aim to extract this transport coefficient from 
experimental heavy-ion data. 

For event-by-event hydrodynamics we show two curves, 
using either the entropy density (blue open circles) or the 
energy density weighted (red solid circles) average eccen- 
tricities to normalize the average final flow (v n ). For the 
cllipticity (Fig. 11) this choice is seen to make a differ- 
ence only in rather central collisions (&<4fm), but for 
the triangularity the differences are significant out to av- 



erage impact parameters probed in minimum bias sam- 
ples, b < 8 fm. As stated earlier, we prefer the energy 
density weighted eccentricities (solid circles) as deforma- 
tion measures because energy density and pressure are 
closely related through the EOS, and it is the pressure 
gradients (and their anisotropics) that drive the collec- 
tive flow (and its anisotropics). 

For the single-shot hydrodynamic simulations, a ques- 
tion arises as to how exactly one should construct the 
ensemble-averaged smooth initial profile which is then 
evolved hydrodynamically. We have explored three rea- 
sonable procedures (variations of which have been used 
in the literature) and show them as dashed lines in 
Figs. 11 and 12. For the lines labeled by stars, we ro- 
tate the entropy density for each fluctuating event by the 
corresponding entropy-weighted participant-plane angle 
V'n P ( s ) ( n ~ 2 5 3, see Eq. (19)), 6 then average the rotated 



Note that for computation of -03 we rotate the events by a differ- 



17 



entropy profiles, compute the eccentricity e n (s) of the re- 
sulting average entropy density profile and convert it to 
energy density using the EOS for input into the hydro- 
dynamic code. For the lines labeled by crosses, we rotate 
the energy density for each fluctuating event (obtained 
from the EOS) by the corresponding energy-weighted 
participant-plane angle xp^ p (e) (see Eqs. (1,2)), compute 
the averaged rotated energy density profile and its ec- 
centricity £ n (e), and use it directly as hydrodynamic in- 
put. For the dashed lines without symbols, finally, the 
averaged initial energy density (and therefore the final 
v n ) are exactly the same as for the lines with crosses, 
but the final v n is scaled by the entropy-weighted (rather 
than energy- weighted) eccentricity of the averaged initial 
profile, where the entropy density is obtained from the 
smooth averaged energy density via the EOS. 

The differences between the different dashed lines il- 
lustrate the uncertainties associated with the choice of 
averaging procedure for the initial state. Keeping in 
mind that a 20% reduction in V2/S2 corresponds (very 
roughly) to an increase of 77/s by l/An [5], one sees that 
these differences are not negligible if one aims for quan- 
titative precision in the extraction of the specific shear 
viscosity. Comparing the three dashed lines we see that 
it doesn't make much difference whether we use the s- 
weightcd or e- weighted participant-plane angles to rotate 
the events before superimposing them (the dashed lines 
without symbols and with stars are all very close to each 
other), but that in the more central collisions we obtain 
significantly different values for the conversion efficiencies 
v n /e n if we normalize by e- or s- weighted mean eccen- 
tricities. Even though they look similar in Fig. 5a, at 
small impact parameters £ pa rt(e) and (£2(e)} arc larger 
than £ P art(s) and (£2(s)}, respectively, and this is the 
main reason why the red and blue lines in Fig. 11 di- 
verge at small b, for both event- by- event (solid lines) and 
single-shot hydrodynamics (dashed lines). 

An apples-to-applcs comparison between event-by- 
event and single-shot hydrodynamics (and between the- 
ory and experimental data) therefore must ensure that 
the same (or at least conceptually compatible) eccentrici- 
ties are used to normalize the anisotropic flow coefficients 
that are to be compared. In Figs. 11, 12 we should there- 
fore compare blue solid with blue dashed, or red solid 
with red dashed lines, but not curves of different colors. 

Even this is not good enough if one wants to accurately 
assess the relative space-to-momentum anisotropy con- 
version efficiency in single-shot and event-by-event hy- 
drodynamics: in the single-shot hydro curves we use £ par t 
to normalize the final elliptic flow, whereas the event-by- 
event hydro results were normalized with (£2) = (£ P art). 
While each of these eccentricity measures makes per- 



ent angle before averaging than for V2, i.e. V3 and V2 are obtained 
from two different single-shot hydrodynamic runs, starting from 
different averaged initial energy density profiles. 



feet sense in its own context, they differ at large im- 
pact parameters, £ par t being larger (see Figs. la,b). To 
avoid this problem we have added in Figs. 11 and 12 
an additional "mixed ratio" (dash-dotted purple line) 
which normalizes the ensemble-averaged anisotropic flow 
(v n ) (n = 2, 3) from event-by-event hydrodynamics (used 
in the ratio (v n )/(e n ) denoted by solid lines with solid 
red circles) by the mean e-wcighted eccentricity e n from 
single-shot hydrodynamics (used in the ratio v n /e n de- 
noted by dashed lines with crosses). This dot-dashed 
purple line agrees almost perfectly with the solid red line 
with circles over most of the impact parameter range, 
except for peripheral collisions with b > 10 fm where e n 
and (e n ) begin to diverge. The red dashed lines with 
crosses and purple dash-dotted lines show the anisotropic 
flows from single-shot and event- by- event hydrodynam- 
ics normalized by the same eccentricity measure charac- 
terizing the fluctuating event sample. Their comparison 
allows an unambiguous assessment of the different effi- 
ciencies of single-shot and event-by-event hydrodynam- 
ics in converting initial eccentricities to final momentum 
anisotropics. Their ratio is shown in Fig. 13. 

From Fig. 13a one concludes that, for ideal hydrody- 
namics, event-by-event fluctuations on average reduce 
the efficiency of the fluid in converting initial source el- 
lipticity into elliptic flow. Over most of the centrality 
range this reduction is about 4% for pions and about 
twice as large for protons, and it is similar for MC- 
KLN and MC-Glauber initial profiles. In very central 
collisions the ratio of conversion efficiencies for event-by- 
event vs. single-shot hydrodynamics is closer to 1, but 
it degrades strongly in very peripheral collisions where 
event- by- event evolution generates on average 30—40% 
less elliptic flow than single-shot hydrodynamics. The 
generic tendency of event- by- event hydrodynamic evolu- 
tion of fluctuating initial profiles to generate less elliptic 
flow than expected from hydrodynamic evolution of the 
corresponding smooth average profile has been observed 
before [27, 30]; our systematic study in Fig. 13a quantifies 
this effect over the full range of collision centralitics. 

The situation with triangular flow, shown in Fig. 13b, 
is quite different: event- by- event propagation of initial- 
state flutuations can lead to an increase or decrease of the 
triangular flow compared to single-shot hydrodynamics, 
depending on particle mass (pions or protons), the nature 
of the fluctuations (MC-Glauber or MC-KLN), and col- 
lision centrality. Contrary to elliptic flow, in peripheral 
collisions event-by-event evolution leads to significantly 
larger average triangular flow than single-shot hydrody- 
namics. 

We expect that non-zero viscosity will dampen fluctua- 
tion effects and somewhat reduce the differences between 
event- by- event and single-shot hydrodynamic evolution 
of elliptic and triangular flow shown in Fig. 13. Nev- 
ertheless it appears that, for quantitative studies of the 
influence of viscosity on the generation of anisotropic col- 
lective flow, event-by-event hydrodynamic evolution is an 
essential and indispensable ingredient. 
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FIG. 13: (Color online) Ratio of the average flow coefficient (v n ) from event-by-event hydrodynamics and the corresponding 
mean v n from single-shot hydrodynamics, as a function of impact parameter in 200 A GeV Au+Au collisions, for n — 2 (a) and 
n~3 (b). Shown are the ratios for directly emitted pions (circles) and protons (diamonds) from fluctuating events using the 
MC-Glauber (dashed lines) and MC-KLN models (solid lines). Average events for computing v n using single-shot hydrodyamics 
were obtained by rotating the energy density of each event by i/^ p (e) before superimposing them. 



C. Elliptic flow fluctuations 

Similar to what is shown in Figs. la,b for the initial 
source ellipticities, Fig. 14 shows the elliptic flow mea- 
sures (v 2 ), ^2{2} and «2{4} from event-by-event hydrody- 
namics, together with v 2 from single-shot hydrodynamic 
evolution of the corresponding averaged initial profile, for 
pions and protons, using MC-Glauber and MC-KLN ini- 
tializations, respectively. f2{2} and ^{4} are defined in 
analogy to Eqs. (9,10) by 

v 2 {2f = (vi), (23) 
« 2 {4} 4 = 2(vl) 2 -(v 4 2 ). (24) 

Here v 2 is calculated event-by-event via Eq. (22) from the 
Cooper-Frye spectrum at freeze-out (with zero statistical 
uncertainties since it is determined with mathematical 
precision by the event-by-event hydrodynamic output). 

As in Fig. 1, open stars show the central values for 
\Z|v2{4} 4 | whenever i>2{4} 4 turns negative, and open- 
ended error bars indicate that the error band for t>2{4} 4 
ranges from positive to negative values. Similar to the 
ellipticities shown in Fig. 1, the latter happens at small 
impact parameters, but for the MC-Glauber model the b- 
range over which this happens for i>2{4} (for both pions 
and protons) is somewhat larger than for e{4}. Still, 
«2{4} is compatible with zero over this entire range, and 
we do not find statistically significant negative values for 
«2{4} at small impact parameters. At large 6>12fm 
v 2 {4} turns negative for both pions and protons when 
we use MC-Glauber initial conditions whereas it remains 
positive for MC-KLN initial profiles. 

By comparing v 2 (open circles in Fig. 14) with (v 2 ) 
(solid green circles) one sees that in mid-central to pe- 
ripheral collisions the ^-suppression from event-by-event 



hydrodynamic evolution is of the same order as or (es- 
pecially for protons) even larger than the difference be- 
tween v 2 {2} and (v 2 ) (solid blue vs. solid green circles) 
that arises from event- by- event flow fluctuations. As a 
result, v 2 {2} from event- by- event hydrodynamics lies in 
peripheral collisions even below v 2 from single-shot hy- 
drodynamics, in spite of its fluctuation-induced enhance- 
ment. 

Similar to Eqs. (11)-(13) we can test whether the v 2 
fluctuations from event to event have Gaussian or Bessel- 
Gaussian distributions. This is done in Fig. 15. The 
upper set of curves (thick lines) test the i>2-analoguc 
of relation (13) whereas the lower set (thin lines) tests 
the validity of Eq. (11). (In the lower set of curves we 
dropped all b- values for which the error band for t>2{4} 
extends ton negative values.) Just as we saw for the ini- 
tial ellipticities in Fig. 3, both the Gaussian and Bessel- 
Gaussian hypotheses for i>2-fluctuations are seen to hold 
quite well in mid-central (4<6<10fm) collisions. The 
Bessel-Gaussian hypothesis breaks down in peripheral 
collisions (b > 10 fm). Whether it holds (as expected [22]) 
in central collisions is a question that, with our present 
statistics of 1000 events per &-value, we can not reliably 
answer, but we see no indications for the opposite. The 
assumption of Gaussian ^-fluctuations breaks down in 
central collisions (6<5fm), as expected. For the MC- 
Glauber model it also breaks down in very peripheral 
collisions, whereas for MC-KLN initial conditions the fi- 
nal elliptic flow exhibits a nice Gaussian distribution all 
the way to the largest impact parameters. 

Overall, a comparison of Figs. 15 and 3 (as well as of 
Figs. 14 and 1) shows that the statistical properties of v 2 
fluctuations are qualitatively similar but quantitatively 
different from those of the initial cllipticity fluctuations. 
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FIG. 14: (Color online) Different measures for the final elliptic flow 112 (similar to Fig. 1) for directly emitted pions (a,b) and 
protons (c,d) as functions of impact parameter from event-by-event ideal fluid dynamics, using MC-Glauber (a,c) and MC-KLN 
(b,d) initial conditions for 200 A GeV Au+Au collisions. 



This is consistent with the fact that the main driver for 
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FIG. 15: (Color online) Similar to Fig. 3, but for the elliptic 
flow V2 of pions and protons. See text for discussion. 



elliptic flow is the initial ellipticity, but that eccentricity 
coefficients of higher harmonic orders affect the evolution 
of V2 weakly but measurably through non-linear mode- 
coupling effects. 

VI. SUMMARY AND CONCLUSIONS 

In this work we presented a comprehensive analysis 
of event- by- event shape fluctuations in the initial state 
and flow fluctuations in the final state of relativistic 
heavy-ion collisions, as quantified by the first four non- 
trivial harmonic eccentricity and flow coefficients, e n and 
v n (n = 2,3,4,5). Using the MC-Glauber and MC-KLN 
models to generate fluctuating initial entropy and energy 
density profiles, we explored the centrality dependence of 
a number of different variants of these anisotropy mea- 
sures that are being used by practitioners in the field, 
and compared them with each other. Although they all 
exhibit similar qualitative behaviour, quantitative differ- 
ences exist and must be carefully taken into account in 
the theoretical analysis of experimental data. As far as 
we know, ours is the first comprehensive analysis quan- 
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tifying these differences for both the Glauber and Color 
Glass Condensate models. 

We list a few key results: 

- The average and mean ellipticities {£2) and £2 agree 
with excellent accuracy over a wide range of impact pa- 
rameters, but diverge in very peripheral collisions (> 60% 
centrality) where £2 > (£2} (both for participant-plane 
and reaction-plane averaged profiles) . 

- The average energy and entropy density weighted ec- 
centricities agree with excellent accuracy over a wide 
range of impact parameters, except for central collisions 
(6<4fm) where (e n (e)) > (e n (s)). 

- Whether the fluctuating entropy density distributions 
for individual events are first converted to energy den- 
sity and then rotated by ^^(e) an d averaged, or first 
rotated by ip pp (s) and averaged and then converted to 
energy density has very little influence on the shape of 
the resulting smooth average initial energy density profile 
for single-shot hydrodynamics. We prefer (and propose 
as standard procedure) the conversion to energy density 
as the first step, since in event- by- event hydrodynam- 
ics the energy density gradients of each event generate 
(through the EOS) the pressure gradients that drive the 
evolution of collective flow. 

- The shortcut of using reaction-plane averaging to gen- 
erate a smooth profile for single-shot hydrodynamics with 
ellipticity approximately equal to £{4} of the ensemble, in 
the hope of generating with a single hydrodynamic run 
an elliptic flow v-i that can be directly compared with 
i>2{4} measurements, works only in the 0—40% central- 
ity range. For peripheral collisions this method cannot 
be trusted. 

- The assumption of Bcsscl-Gaussian fluctuations for ini- 
tial source ellipticity and final elliptic flow work well for 
b< lOfm but breaks down in more peripheral collisions. 
For more peripheral collisions the hypothesis that £2 and 
i>2 are Gaussian distributed works better than the Bessel- 
Gaussian assumption, but it breaks down for b < 5 fm. 
For MC-Glaubcr initial conditions, directly emitted pi- 
ons and protons feature negative values of i>2{4} 4 in very 
peripheral collisions. The fluctuations of initial source el- 
lipticities and final elliptic flow values have qualitatively 
similar but quantitatively different statistical properties. 

- Except for rather central collisions, the eccentricities 
(£2), (£4) and (£5) from the MC-KLN model are all sig- 
nificantly larger than those from the MC-Glauber model. 
In contrast, (£3) is numerically very similar for the two 
models over most of the impact parameter range. The 
viscous suppression of triangular flow V3 may thus allow 
for a determination of the QGP shear viscosity (?//s)qgp 
that is free from the large model uncertainties that arise 
from the different MC-Glauber and MC-KLN ellipticities 
when using vi for such an extraction [49]. 

- The second and fourth order eccentricities £2 and £4 
are strongly correlated by collision geometry, and D4 re- 
ceives strong contributions even from a purely elliptical 
deformation of the final flow velocity distribution. These 



complications make 1)4 a poor candidate for systematic 
studies of viscous effects on the evolution of collective 
flow. Similar comments apply to W5 since it couples via 
mode-coupling to triangularity from fluctuations and to 
ellipticity from collision geometry. This mixture of con- 
tributions from conceptually different origins complicates 
a systematic analysis. In general, flow coefficients v n of 
high harmonic order (n > 3) show poor correlation with 
the eccentricity coefficients e n of the same harmonic or- 
der, except for very central collisions where all eccentric- 
ities are driven by fluctuations alone (and not by overlap 
geometry) . 

- In spite of non- linear mode-coupling effects, the basic 
response of elliptic flow v 2 to ellipticity £2 , and of trian- 
gular flow ^3 to triangularity £3, is approximately linear. 
These two observables thus remain prime candidates for 
systematic studies of viscous effects on collective hydro- 
dynamic flow. 

- Event-by-event hydrodynamics generates harder px- 
spectra for the emitted hadrons than single-shot hydro- 
dynamic evolution of the corresponding averaged initial 
profile. This is due to additional radial flow generated by 
large pressure gradients arising from "hot spots" in the 
initial fluctuating density distribution. The hardening 
effect is particularly strong in peripheral collisions which 
produce small fireballs that fluctuate strongly. 

- Evcnt-by-event hydrodynamic evolution of fluctuating 
initial conditions leads to smaller average elliptic flow 
than obtained by evolving the corresponding averaged 
initial condition in a single shot. This suppression de- 
pends somewhat on collision centrality, and for ideal flu- 
ids it is generically of order 4-5% for pions and 8-10% for 
protons. The effect is sufficiently large to lead to a sig- 
nificant over-estimate of the fluid's specific shear viscos- 
ity if one extracts it from elliptic flow measurements by 
comparing with single-shot hydrodynamic simulations. 
Even though we expect the discrepancy between event- 
by-event and single-shot hydrodynamics to decrease a bit 
in viscous fluid dynamics, we believe that a quantitatively 
precise experimental determination of i]/s from collective 
flow data will require comparison with cvcnt-by-cvcnt hy- 
drodynamical calculations. 
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FIG. 16: (Color online) Correlation between e„ and e' n , for n = 2, 3, 4 (panels (a-c)). The blue dots are bin averages for bins 
that contain more than 10 events. The thick black lines are linear fits. 
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FIG. 17: (Color online) Event-by-event correlation between the participant plane angles associated with r 2 - and r n -weighted 
eccentricities for harmonic orders n — 3, 4, and 5 (panels (a)-(c)). 



Appendix A: Comparison between eccentricities 
defined with r 2 and r n weights 

We here present a brief comparison between the r 2 - 
weighted eccentricity coefficients e„ (Eq. (17)) and the 
r n -weighted e' n (Eq. (17)), as well as their associated an- 
gles ifj FF and ipn P ■ Fig- 16 shows a scatter plot of e' n vs. 
e n for n = 3, 4, 5. One observes approximate proportion- 
ality (£3 w 1.22 £3, £4 1.48 £4, £5 « 1.80 £5) over most of 
the eccentricity range, with slopes that increase with n. 
So where Fig. 5 shows a decrease of £„ with increasing n 
at large impact parameters, the same is not true for the 



e' n [20] . On the other hand, the linear relations between 
e' n vs. £„ imply that the relations between v n and e' n will 
look qualitatively the same as those between v n and e n 
in Fig. 9, with appropriately rescaled horizontal axes. 

At the same time the participant plane angles associ- 
ated with r 2 -weighted and r™-weighted eccentricities are 
tightly correlated, as shown in Fig. 17. For given n, the 
angles ip FF and ip FF fluctuate around each other, with 
a variance that increases with n, on account of the de- 
creasing values of e n . From a practical point of view, 
we therefore consider both definitions as equivalent, and 
choosing between them is a matter of personal preference. 
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